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ABSTRACT 

We have observed the pulsar in the Crab Nebula at high radio frequencies and high time resolution. 
We present continuously sampled data at 640-ns time resolution, and individual bright pulses recorded 
at down to 0.25-ns time resolution. Combining our new data with previous data from our group and 
from the literature shows the dramatic changes in the pulsar’s radio emission between low and high 
radio frequencies. Below about 5 GHz the mean profile is dominated by the bright Main Pulse 
and Low-Frequency Interpulse. Everything changes, however, above about 5 GHz; the Main Pulse 
disappears, the mean profile of the Crab pulsar is dominated by the High-Frequency Interpulse (which 
is quite different from its low-frequency counterpart) and the two High-Frequency Components. We 
present detailed observational characteristics of these different components which future models of the 
pulsar’s magnetosphere must explain. 

Subject headings: pulsars: general; pulsars: individual (Crab pulsar) 


1. INTRODUCTION 

Instantaneous pulsar radio emission is strongly vari¬ 
able. Its intensity can fluctuate on timescales as short 
as microseconds or even nanoseconds. The mean pro¬ 
file, however, is quite stable, and “components” of the 
mean profile are easy to identify. That is, strong radio 
pulses occur only at certain rotation phases, so that their 
sum produces a clear component in the time-averaged 
mean profile. These components therefore represent lo¬ 
calized, long-lived “hot spots” in the magnetosphere, 
where plasma conditions are favorable for the creation of 
coherent radio emission. The fact that the high-energy 
mean profiles also show distinct components suggests 
that similar active regions within the magnetosphere cre¬ 
ate pulses at high photon energies. 

The major question is, where in the magnetosphere are 
these hot spots? Their locations depend critically on the 
physics of the magnetosphere: unshielded electric fields, 
plasma streams, pair creation, plasma and/or MHD in¬ 
stabilities can all contribute to bursts of radio and high- 
energy emission. Are the hot spots restricted to low al¬ 
titudes, in the open field line regions above the magnetic 
polar caps? Do they extend for large distances, out to or 
even beyond the light cylinder, along high-altitude caus¬ 
tics? Do they exist elsewhere in the magnetosphere, in 
regions not predicted by any current model? 

In this paper we focus on radio emission from the pul¬ 
sar in the Crab Nebula. Although this star has been 
widely studied at wavelengths from low radio frequencies 
through optical, X-rays and y-rays, no extant model of 
that pulsar’s magnetosphere can address all of the avail¬ 
able data. In this paper we present new observations at 
high radio frequencies, and use these to characterize the 
pulsed radio emission from the Crab pulsar in as much 
detail as possible. We present these results as metrics 
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against which present and future magnetospheric mod¬ 
els should be tested. 


1.1. Two Strong Emission Components 

At low radio frequencies (below ~ 5 GHz) and at high 
energies (optical, X-ray, y-ray), the mean profile of the 
Crab pulsar shows two strong features, generally known 
as the Main Pulse and the Interpulse. The phase sepa- 



the same phase strongly suggests the regions emitting 
both components are located close to each other in the 
star’s magnetosphere. 


1 . 1 . 1 . Radio Models 

In models of pulsar radio emission, the Main Pulse and 
Interpulse are usually interpreted as low-altitude emis¬ 
sion from open field lines above the two magnetic poles 
of an orthogonal rotator, seen at a viewing angle ~ 90° 
between the r otatio n axis a nd sight line. Although |Wel-| 
tevrede & Johnston (2008) noted there are too few puL 
sars with interpulscs in the observed population to be 
consistent with a random distribution of viewing angles, 
the orthogonal-rotator model is supported by polariza¬ 
tion obse rvations of some p ulsars with interpulses. For 
instance, Keith et al. (2010) show that the swing of po¬ 
larization position angle across the Main Pulse and In¬ 
terpulse in five pulsars with interpulses is consistent with 
the rotating-vector model (Radhakrisnan & Co oke|196 9) 
and a nearly orthogonal viewing angle. However, this 
model does not easily fit the Crab pulsar, in which the 
polarization position angle shows no variation with p ulse 


phase across the Main Pulse or the Inter pulse (Moffett 
& Hankins 1996[ Slowikowska et al. 2014). In addition 


the rotating-vector model is further called into question 
for the Crab pulsar by the shape of the X-ray torus 
within the Cra b Ne bula, which requires a viewing angle 
~ 60° (Ng & Romani 2004). Because this angle is too 
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far from orthogonal to allow the low-altitude open field 
lines from both magnetic poles to pass within our sight 
line as star rotates, the model of low-altitude polar-cap 
emission seems not to apply to this pulsar. 

1 . 1 . 2 . High-Energy Models 

Models of high-energy pulsar light curves paint quite 
a different picture. Pulsar X-ray and 7 -ray light curves 
often have two well-defined, sharp peaks s eparat ed by 
130° to 180° in pulse rotation phase (e.g., Abdo et al. 


2013). Such light curves are too common to be consis- 


tent with low-altitude emission from the magnetic poles 
of orthogonal rotators. In addition, the highest-energy 7 - 
rays must be emitted from well above the star’s surface, 
because the large pair-creation opacity close to star’s sur¬ 
face precludes t heir escape from low-a ltitude emission re¬ 


gions (e.g., Romani & Watters 2010). Models therefore 
place the emission sites at high altitudes, typically a sig¬ 
nificant fraction of the distance to the light cylinder, in 
an outer gap or an extended slot gap. Field line sweep- 
back and photon pile-up on caustics then creates two 
broad pulses (main pulse and interpulse) when seen at a 
large range of viewing and inclination an gles (e.g., Ro- 
mani & Watters|2010 Harding et al.|2011 ). For the Crab 
pulsar, the phase coincidence of the Main Pulse and In¬ 
terpulse in both radio and high energy bands suggests the 
radio and high-energy emitting regions are regions spa¬ 
tially close in star’s magnetosphere. This also seems to 
disagree with the low-altitude polar cap emission model 
of radio emission. 

1 . 1 . 3 . Alternatives and Complications 

Of course, alternatives to the two -pole, orthogo nal ro¬ 
tator model have been suggested. [Manchest er fe Lyne| 
(1977) argued that that an interpulse is just an ex¬ 
treme version o f a two-peaked pul se aris ing from one 


magnetic pole. Hankins & Cordes (1981) pointed out 
difficulties with the two-pole model of interpulses for 
PSR B0950+08, and noted that both one-pole and two- 
pole models are possib le for that star. The annular gap 
model (e.g., Du et al. 2012) has both radio and high- 
energy emission "coming from a mid-altitude region above 
only one pole. Similarly, [Ravi, Manchester & Hobbs] 
( 2010 ) suggested that both radio and high-energy emis- 


sion m young, fast-spindown pulsars, comes from fan 
beams emitted at high altitudes in the magnetosphere. 
Petrova (2009) proposed that induced Compton scatter¬ 
ing of Main Pulse emission creates all other components 
in the mean profile of the Crab pulsar. 

Furthermore, the Interpulse of the Crab pulsar changes 
dramatically at high radio frequencies. Abov e ~ 5 GH^ 
it shifts to ~ 7° earlier in rotation phase (Moffett & 


Hankins|1996 ). The spectral and temporal nature of sin¬ 


gle Interpuiscs is very different from that of single Main 


Pulses (Hankins & Eilek 


2007). We consider the Low- 


Frequency Interpulse and High-Frequency Interpulse to 
be two totally separate components. Even though the 
two Interpulses are close to each other in phase, and both 
fall c omfort ably within the broad high energy Interpulse 
(e.g.,|Abdo et al.][ 2010 |), the disparity between the radio 
characteristics of the two components suggests they arise 
from different regions within the pulsar’s magnetosphere. 

1 . 2 . Even More Radio Components 


TABLE 1 

Components of the Mean Profile 


Component 

Acronym 

Frequency 

Range 

Precursor 

PC 

0.3 - 

■ 0.6 

GHz 

Main Pulse 

MP 

0.3 - 

4.9 

GHz 

High-Frequency lnterpulse 

HF IP 

4.2 - 

■ 28.4 

GHz 

Low-Frequency lnterpulse 

LF IP 

0.3 - 

■ 3.5 

GHz 

High-Frequency Component 1 

HFC1 

1.4 - 

■ 28.0 

GHz 

High-Frequency Component 2 

HFC2 

1.4 - 

■ 28.0 

GHz 

Low-Frequency Component 

LFC 

0.6 - 

■ 4.2 

GHz 


Note. — Frequency ranee over which compon ent is de¬ 
tected in mean profiles l|Moffett. & Hankins 1996] and this 
work). Occasional single pulses may be detected outside this 
range, but they are too rare to contribute to the mean profile. 
These acronyms are used only in the Tables and Figures. 


This is not yet the entire story of radio emission from 
the Crab pulsar. Its mean radio profile is complex, with 
seven separate components identified so far. We list these 
in order of incr easin g p hase in Tab le [Tj following the 
nomenclature of Moffett & Hankins (1996). 

The Crab pulsar is unusual in having such a large num¬ 
ber of components spread throughout the rotation pe¬ 
riod. While there are many pulsars whose mean profiles 
contain multiple co mponents (e.g., PSR B123 7+25 which 
has five components, Hankins & Wright|1980 1 those com¬ 
ponents are typically bunched closely together in phase 
and are generally considered part of a single mean-profile 
component (for instance within the Main Pulse). By 
contrast, the seven components found in the Crab pul¬ 
sar are widely spaced in phase, extending throughout the 
star’s full rotation period. There are also a few pulsars, 
such as PSR B0826—34 and PSR B1929+10, which emit 
over most or all of their rotat ion period. While some 
authors (e.g., Biggs et al. 19851 have suggested the mag¬ 
netic axis of such stars is oriented close to the line of sight 
(so-called aligned rotators), this picture is not necessar¬ 
ily co nsistent with th e observed p osition angle behavior 
(e.g., Rankin & Rathnasree 1997). 

While the mapping from magnetospheric position to 
rotation phase is far from linear, the wide phase spread 
of the Crab’s components suggests that many separate 
radio emission sites exist throughout the star’s magneto¬ 
sphere. We have not seen much discussion in the liter¬ 
ature of where in the magnetosphere these components 
may arise. It has occasionally been suggested that the 
Precursor is generated at low altitudes near the pul¬ 
sar su rface in the open- field-line region of the polar cap 
(e.g., Lyne et al. 2013). This suggestion agrees with 
caustic models which locate Main Pulse emissi on, both 


radio and high-energy, at high altitudes (e.g., Harding 


et ah 2011 ); low-altitude emission arising close to fEe 


polar cap will lead the high-altitude Main Pulse emis¬ 
sion b y sev eral degrees of ph ase. Both Hankins fe Eilek| 
(2007]) and lLyutikov| (|2008j) speculated that the High- 
Frequency lnterpulse might arise from high altitudes, 
close to the light cylinder, but these ideas have not been 
developed far enough to test against mean profiles. We 
are not aware of any suggestions for the spatial locations 
of the other components (the two High-Frequency Com¬ 
ponents, and the Low-Frequency Component). While 


recent studies of pulsar magnetospheres (e.g., Li et al. 
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(12012 ) ; Kalapotha rako s 
Kalapotharakos 2Ul(jb ) 


et al. (2012); Contopoulos & 


alapotharakos (2010)) do not specifically model the ra- 
dio prohle, the extended current sheets and dissipation 
regions revealed by their simulations seem very likely to 
be sources of coherent radio emission. 


1.3. Our Work: Mean Profiles and Single Pulses 

Our primary focus in this paper is the phenomenology 
of the mean radio profile of the Crab pulsar. We present 
continuously sampled data, recorded between 9 and 43 
GHz at sub-microsecond time resolution. We use these 
data to derive and characterize mean profiles at these 
high radio frequencies, and also identify and character¬ 
ize the statistics of bright single pulses within the data 
streams. In this paper we use the term “mean profile” 
to describe the intensity as a function of rotation phase. 
This quantity is derived by summing the received inten¬ 
sity synchronously with the star’s rotation period. The 
same data product is also often called the “light curve”, 
especially in the high-energy community. 

We have also recorded strong individual pulses between 
0.33 and 43 GHz, at time resolution down to a frac¬ 


tion of a nanosecond (Hankins & Eilek 


et al. 20101. These occasional strong pulses emitted by 


2007 Crossley 


the Crab p ulsar led to its discovery (Staelin & Reifen- 


stem 


1968) lon g before its periodicity was determined 


\ Cornelia et af. 1969 ). These pulses have been called “gi¬ 


ant pulses”, because their intensities are far greater than 
the typical pulse. However, there is no standard cri¬ 
terion for labeling these single pulses as “giant”, and in 
fact the distribution of pulse amplitudes in the Crab pul- 
sar appears to be continu ous from weak to strong, (e.g., 


Karuppusamy et al. 2010). In this work we will there¬ 
fore",] ust refer to “single pulses”, meaning those bright 
enough to be clearly discernible above the system noise. 
Such pulses are relatively rare; in most rotation periods 
the pulses are below the noise level. When we do de¬ 
tect strong, single pulses, they occur at the phase of the 
Main Pulse and the Interpulse, and very occasionally at 
the phases of the High-Frequency Components. 

In this paper we present the ensemble characteristics 
of the Crab pulsar’s radio emission revealed by our data. 
We describe our observations in Section [2} we extend 
up to 43 GHz previo us s tud ies at l o wer r adio frequen¬ 
cies, such as that of Cordes et al. (2004|). In Section 


f we present the composite results: mean profiles as a 
metion of frequency, and the phases and widths of the 
components of the mean profile. In Section [4] we present 
ensemble characteristics of individual Mean Pulses and 
Interpulses: their arrival phase and width, their fluence 
distributions, and the “burstiness” of arrival times. In 
Section [5] we summarize our results and discuss the con¬ 
straints they place on current and future models of the 


In a separate paper to follow, (Hankins et al. 2015 

Pa- 

per 2) we extend our previous sing 

le-pulse work ( 

Hankins 

et al. 2003 

Hankins & Eilek|2007 

Crossley et al. 

2010 

) to 


new observations of bright single pulses at the phases of 
the Main Pulse, the Low-Frequency and High-Frequency 
Interpulses, and the two High-Frequency Components. 


2. OBSERVATIONS 


We recorded pulses from the Crabpulsar in observing 
sessions at the Arecibo Observatorjr] between 2002 and 
2009, the Robert C. Byrd Green Bank Telescope (GBT)(^] 
in 2009, 2010, and 2011, the Jansky Very Large Array 
(VLA) from 1993 to 1999. 

We used two separate data acquisition systems. To ob¬ 
tain continuously sampled data (presented in Sections [3] 
and [4| , we used the Green Bank Ultimate Pulsar Pro¬ 
cessing Instrument (GUPPI) at the GBT, using receiver 
bands 8-10, 12-15.6, 18-22.4, 22-26.5 and 41-46 GHz and 
the single linear polarization available at 28-32 GHz. We 
operated the GUPPI system in the “coherent search” 
mode, in which the 800-MHz bandwidth was split into 
32 frequency bands. Each 25-MHz band was then coher¬ 
ently dedispersed at the nominal Jodrell Bank dispersion 
measure^ Total intensity samples were then formed from 
each time series and every 16 samples were accumulated 
before being written to disk, resulting in a final time res¬ 
olution of 640 ns. 

To capture individual puls es at higher time resolution 
(as reported in Section |4.3| and in more detail in Pa¬ 


per 2), we used our Ultra High Time Resolution System 
(UHTRS) at the Arecibo telescope, the GBT and the 
VLA. In this system, a large-memory digital oscilloscope 
sampled and recorded the received voltages of both polar¬ 
izations. A square-law detector with 20-/zs time constant 
was used to establish a signal threshold, typically set to 
6 times the smoothed off-pulse root-mean-square noise 
level. When the total intensity exceeded the threshold, 
the oscilloscope was triggered to sample and digitize the 
voltages at appropriate Nyquist rates up to 10 Gigasam- 
ples per second during a sampling window centered on 
the pulse component of interest. For the observations 
at the Arecibo Observatory in 2009 we also used a re¬ 
altime digital dedisperser based on the CASPEB^ iBob 
device as the oscilloscope trigger generator. The 8-bit 
oscilloscope-sampled data were transferred to comp uter 
disk for subseq uent off -lin e co h eren t dedispersion |Han-| 


kins 


1971 Hankins & Ric kett| [l975), which allowed time 


resolution up to the inverse of the receiver bandwidth, 
about 0.2-0.4 ns. During the data transfer time from 
oscilloscope to computer disk, 10-30 s, the data acquisi¬ 
tion wa s di sabled, and no pulses could be captured. In 
Section l4~3l we use data from the UHTRS to characterize 
widths of single pulses; in Paper 2 we present much more 
single pulse data taken with the UHTRS. 

3. CHARACTERISTICS OF THE MEAN PROFILE 
3.1. Mean Profiles as a Function of Frequency 
We have exten ded the radio mean profiles of [Moffett 


& Hankins (1996) up to 28 GHz in order to test for the 


frequency dependence of the spacing and width of the 


5 The Arecibo Observatory is operated by SRI International un- 
der a cooperative agreement with the National Science Foundation 
(AST-1100968), and in alliance with Ana G. Mendez-Universidad 
Metropolitana, and the Universities Space Research Association. 

6 The Robert C. Byrd Green Bank Telescope and the Jansky 
Very Large Array (JVLA) are instruments of the National Radio 
Astronomy Observatory, facilities of the National Science Founda¬ 
tion operated under cooperative agreement by Associated Univer¬ 
sities, Inc. 

7 The Jodrell Bank Crab Pulsar Monthly Ephemeris, 
http: / / www.jb.man.ac.uk/pulsar/crab.html 

0 Collaboration for Astronomy Signal Processing and Electronics 
Research, https://casper.berkeley.edu 
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Pulse Phase (Degrees) 


Fig. 1.— Mean profiles for a number of frequencies are shown with the formal Gaussian fits to relevant component s over-plotted . The 
Low-F requency Component, LFC , and the High-Frequency Components, HFC1 and HFC2, were first identified by Moffett & Hankins 
fl996|. The Main Pulse, MP , Low-Frequency Interpulse, LFIP, and High-Frequency Interpulse, HFIP , are also labelled. The O.T3-GHz 
and 0161-GHz profiles show the Precursor, PC, which only appears below about 1 GHz. The pulse phase of the peaks of the fitted Gaussians 
are denoted by vertical b ars. The time resolut ion for each pro file is shown by th e small horiz ontal bar at the left of each profile. Sources: 
0.43, 1.18, 3.5, 4.1 5 GHz: |Cordes et a 1. | :]2001 [i : 0.61, 1.4 GHz: |Lyne et al.| (20131; 0.812 GHz: |Lundgren et al.| l| 1 99o[l : 4.885 GHz: |Moffett| 
|&; Hankins] Q1999} ; 8.9 GHz: S. Hansom, private communication CZU 11 ): u > 10 GHz: this work - 


mean profile components. In Figure [l] we show our new 
profiles aligned with a selection of previously published 
profiles. As insufficient timing information was available, 
and because the observations were made at widely sep¬ 
arated epochs, we have made two choices to align the 
profiles. For frequencies <, 5 GHz, we have set the Main 
Pulse phase at 0°. For frequencies ,> 5 GHz we have as¬ 
sumed that the High-Frequency Interpulse phase is fixed 
at 137.56° after the center of the Main Pulse. This phase 
separation is based on several data sets which show a de¬ 
tectable Main Pulse as well as a High-Frequency Inter¬ 
pulse at v > 5 GHz. Although the mean profile of the 
Main Pulse is very weak above ~ 5 GHz, we have de¬ 
tected occasional single pulses at the phase of the Main 
Pulse a t all frequencies observed, up to 43 GHz (as in 
Section |4Tj). 

Figure [Tj shows that the nature of the mean profile 
changes dramatically between low and high radio fre¬ 
quencies, as follows. 


• The mean profile shape changes abruptly at ~ 5 
GHz. The Main Pulse, which dominates the mean 
profile at low radio frequencies, nearly disappears 
above ~ 5 GHz. The high-frequency mean profile 
is dominated by the Interpulse and the two High- 
Frequency Components. 

• The Interpulse undergoes a sudden, discontinuous 
phase shift of ~ 7° around 5 GHz. This transition 
coincides with the frequency above which the spec¬ 
tral and temporal characteristics of the Interpulse 
change dramatically (Hankins & Eilek 2007). 

• The two High-Frequency Components can be iden¬ 
tified in the mean profile from 1.4 GHz to 28 GHz. 
Both components drift in phase, appearing at later 
phases at higher frequencies. 

• The Precursor can be identified in mean profiles at 
0.43 and 0.61 GHz, but is not seen at any higher 
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frequency. Similarly, the Low-Frequency Compo¬ 
nent can be seen in mean profiles from 0.61 to 4.2 
GHz, but is not seen at higher frequencies. 


3.2. Phases of Components in the Mean Profile 

We simultaneously fitted single Gaussians to each of 
the profile components in Figure |T1 ignoring components 
that are below the profile noise level. The phase of a 
component is defined as the peak of the fitted Gaussian. 
The phase uncertainties are the quadrature sums of the 
ler Gaussian centroid p ositions as obtaine d from the for¬ 
mal fitting procedure (Press et al. 1986). The sample 
point weighting in our least-squares fitting procedure is 
based on the signal-to-noise ratio, so that profiles with 
higher S/N ratios dominate the fits. 



Frequency (GHz) 

Fig. 2.— The phases of the fitted components of the mean pro¬ 
files. The lcr uncertainty error bars for the formal phase fits are not 
shown when they are smaller than the plot symbol size. The lines 
are least-squares fits of the form (j> = a + 61og(z/), where </> is the 
pulsar rotational phase in degrees, and v is the frequency in GHz. 
The fit coefficients are given in Table [2] The dotted lines indicate 
that the Main Pulse phase is defined tobe zero for 0.4 < v < 5 GHz 
and the High-Frequency Interpulse phase is defined to be 137.56° 
for 5 < v < 28 GHz. We did not include the 0.61-GHz Low- 
Frequency Component in the formal fit as its position is biased by 
the stronger, overlapping Precursor. 


In Figure [2] we show the p hase s of each of the fitted 
components. As in Section |3.1[ we assume the Main 
Pulse has a phase of <^>mp = 0 , and the High-Frequency 
Interpulse has a phase of 0 hfip = 137.56°. For the two 
High-Frequency Components we computed a fit of phase 
vs. log(frequency), = a + b logiz, with v measured in 
GHz. We tried fitting polynomials of the form (j> = Co + 
C 1 P + C 2 JZ 2 but found that the y 2 error of the polynomial 
fit to be much larger. We present the fitting coefficients 
a and b in Table [2] along with the extrapolated values 
of the phases at a reference frequency of 4 GHz. Before 
fitting for ^lfip(^) and (6 lfc(p) we extrapolated their 
observed epoch phases to MJD 530 00 usin g freq uency 
interpolated values of dfi/dt from Lyne et al. (2013). The 
resulting secular phase changes are" in ail cases smaller 
than our fitting uncertainties. 

These results quantify the trends apparent in Fig¬ 
ure [T] The Low-Frequency Interpulse shows very little 
phase drift with frequency. The phase separation be¬ 
tween the Low-Frequency Interpulse and the Main Pulse, 


TABLE 2 

Mean Profile Component Phases 


Component 

a 

(degrees) 

b 

(degrees) 

0(4 GHz) 
(degrees) 

MP 

0.0 

0.0 

0.0 

HF IP 

137.56 

0.0 

137.56 

LF IP 

145.36 ± 0.03 

-0.77 ± 0.12 

144.9 ±0.2 

HFC1 

192.9 ±0.8 

24.2 ±0.8 

207.5 ± 2.5 

HFC2 

255.7 ±0.4 

15.3 ±0.4 

264.9 ± 1.4 

LFC 

322.46 ±0.09 

1.7 ±0.5 

323.1 ± 0.8 


Note. — The component phases are fitted by the ex¬ 
pression <t>(y) = a + b logiq with frequency u in GHz. The 
uncertainties in a and b are lcr; uncertainties in o(4 GHz) 
are from quadrature sums. The phases of the Low-Frequency 
Component and the Low-Frequency Interpulse have been ex¬ 
trapolated to MJD 53000. 


measured from our mean profiles, varies by only ~ 1.1° 
between 430 MHz and 3.5 GHz. Similarly, inspection of 
our single-pulse data showed no evidence for phase drift 
with frequency between the Main Pulse and the High- 
Frequency Interpulse. Except for the ~ 7° phase jump 
when the Low-Frequency Interpulse disappears and the 
High-Frequency Interpulse appears, each Interpulse re¬ 
mains at a nearly constant phase relative to the Main 
Pulse. 


Hankins & Fowler ( 1986]) measured the separation be¬ 
tween the Main Pulse and the Interpulse from 0.196 to 
2.695 GHz. Their phase data and frequency-dependent 
fits to the separations are consistent with the current 
work. They obtained ^lfip(1 GHz) = 145.9 ± 0.6° with 
a frequency exponent of b = —0.68 ± 0.40. This com¬ 
pares well with our </>lfip(1 GHz) = 145.36 ± 0.06° with 
a frequency exponent of b = —0.77 ± 0.12. The secu¬ 
lar change of </>lfip found by |Lyne et al.| (|2013|) over the 
approximately 20 years from Hankins and Fbwler’s obser¬ 
vations to Lyne’s MJD 53000 reference date is far smaller 


than the quoted separation estimation error of Hankins 


& Fowler (1986). 

In the frequency transition region between the Low- 
Frequency Interpulse and the High-Frequency Interpulse, 
we measured at 4.15 GHz, the Main Pulse to Interpulse 
separation as 139.1° ± 0.5°, intermediate between the 
phases of the Low-Frequency and High-Frequency Inter¬ 
pulses. Inspection of single-pulse data in this frequency 
range shows that both types of Interpulses can occur 
around ~ 4 GHz; we thus interpret the intermediate 
phase of the mean-profile Interpulse at 4.15 GHz as the 
result of this mixture. 

Similarly, neither the Precursor nor the Low-Frequency 
Component show any strong phase drift over the fre¬ 
quency range in which they appear in our mean profiles. 
However, the phases of the two High-Frequency Com¬ 
ponents increase with frequency; at 28 GHz both High- 
Frequency Components lag their 1.4-GHz counterparts 
by more than 20°. 

The approximate constancy of the Main Pulse- 
Interpulse phase sepa ration agrees with sever al other in¬ 
terpulse pulsars (e.g., 
those authors found t 


Hankins & Fowler) 19861. However, 
lat the Main Puise-lnterpulse sepa¬ 
ration in PSR B1944+17 changes significantly, by ~ 15° 
between 0.43 and 2.4 GHz, reminiscent of the phase drift 
we find in the High-Frequency Components of the Crab 
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pulsar. Based on the phase shift in PSR B1944+17, 
in addition to polarization anomalies and syn chronous 
nulling in th e Main Pulse and Interpulse, Hanki ns fe 
Fowler (19861 speculated the radio emission is either from 
aTsingle magnetic pole or from high altitudes in that star. 
Might this also be the case for the Crab pulsar? 


3.3. Widths of Components in Mean Profile 

The Gaussians we fitted to the components in the mean 
profile for phase determination can also be used to char¬ 
acterize the width of each component. We show in Fig¬ 
ure [3] the measured full-width at half-maximum of each 
fitted Gaussian, as a function of frequency. To quan¬ 
tify the relation between width, w, and frequency, we 
fitted a linear function in log space to the data points: 
log w(v) = a+ft log v, where v is the observing frequency 
in GHz. We present the fitting coefficients a and j3 in 
Table |3l along with the widths at 4 GHz, evaluated as 
w (4 GHz) = 10 a vP. 

The measured widths of the components are, of course, 
affected by data acquisition systematics, which vary from 
observation to observation. These include the instru¬ 
mental time resolution, the dedispersion technique and 
the value of dispersion measure used, the accuracy of 
the pulsar ephemeris used to calculate the pulse periods 
used to form the mean profiles, and, particularly below 
1 GHz, the variable scattering broadening in the Crab 
Nebula and the interstellar medium. An upper estimate 
of the component width broadened by scattering can be 
obtained by scaling th e values for the scat tering decay 
times, td, measured by Kuzmin et al. (2008) at 111 MHz 
over 2.5 years, as iz -4 , and multiplying the result by 2.46 
to convert td to w, the component full-width at half¬ 
maximum. (Although some authors assume Kolmogorov 
turbulence for the general ISM, which would imply 
the low-frequency pulse broadening of the Crab pulsar is 
m ore consisten t with Gaussian turbulence; e. g., Cross- 
ley et ah[2010). Using the largest value of their observed 
range, 8 < td < 26 ms, we get for our lowest observation 
frequency w(430 MHz) « 2.46 * 26 * (430/111)" * * * 4 = 0.28 
ms, or about 3.1° of pulse phase, which is slightly larger 
than the size of the plot symbols in Figure [3] Since none 
of our observations were made during a known scattering 
event, we conclude that our width measurements are not 
strongly biased by scattering broadening. 

Despite these uncertainties, we find there is very lit¬ 
tle change of the component widths over the observed 
frequency range. The /3 coefficients we find are signif¬ 
icantly non-zero only for the Low-Frequency Interpulse 
and the High-Frequency Component 2, and we note that 
the formal fit to the Low-Frequency Interpulse is strongly 
influenced by two very high S/N points at 0.43 and 1.4 
GHz. Thus, only High-Frequency Component 2 shows 
clear evidence for width changing with frequency. 

For many pulsars the mean profile componen ts are 
broader at low frequencies. For instance, Rankin (19831 
fitted the half-power component widths below ~ 1 GHz 
in several pulsars as w(y ) oc , with /3 ranging from —0.5 
to nearly zero. This “radius-to-frequency mapping” is 
often interpreted as lower-frequency radio emission com¬ 
ing from higher altitudes within the open field line region 
(but still close to the star’s surface). At higher frequen¬ 
cies, however, component widths tend to stabilize (e.g., 


TABLE 3 

Mean Profile Component Widths 


Component 

10 a 

(degrees) 

0 

w (4 GHz) 
(degrees) 

MP 

4.0 ± 1.1 

-0.05 ±0.10 

3.7 ± 1.0 

HF IP 

7.5 ± 1.5 

0.06 ±0.15 

8.2 ±3.1 

LF IP 

4.9 ± 1.0 

0.21 ± 0.03 

6.6 ± 1.3 

HFC1 

25.9 ± 5.0 

-0.03 ±0.08 

24.8 ± 5.0 

HFC2 

10.5 ± 1.1 

0.19 ±0.03 

13.5 ± 1.2 

LFC 

8.1 ± 1.1 

-0.23 ±0.19 

5.9 ± 1.9 


Note. — The component widths are fitted by the ex¬ 
pression log w(u) = a-\-(3 log i/, with v in GHz. The un- 
certanties in a and /3 are lcr; uncertainties in w (4 GHz) 
are from quadrature sum. 


Rankin|[l983 Thorsett~||1991 ). Our finding of frequency- 


independent widths tor most mean-profile components of 
the Crab pulsar agree with this general trend. We are not 

aware of any other pulsar with a mean-profile component 
that broadens at higher frequencies, as High-Frequency 
Component 2 does in the Crab. 
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Fig. 3.— The full-width at half-maximum of Gaussian fit 
widths of the mean profiles in Figure |T] for the Main Pulse, the 
Low-Frequency and High-Frequency Interpulses, the two High- 
Frequency Components, and the Low-Frequency Component. The 
fit components are given in Table [3] The error bars are not plotted 
when the fitted width uncertainty is less than 10% of the fitted 
width. For comparison with single-pulse widths (in Figurefflj). note 
that 1° of phase ~ 90 fis for the Crab pulsar. 


4. ENSEMBLE CHARACTERISTICS OF SINGLE PULSES 
4.1. Single Pulse Arrival Phases 

In addition to the mean intensity profile, we used the 
GUPPI system to track the occurrence rate of pulses 
strong enough to be detected as single pulses (above the 
system noise, typically 6cr). In Figured w e com bine our 
new d ata (above 5 GHz) with data from [Cordes et ah 
(20041 to show the frequency evolution of pulse-count his- 


tograms of the arrival phase of individual bright pulses. 
This data product highlights strong, infrequent single 
pulses which may be hidden in, or distributed differently 
from, the mean intensity profile (Figure [Tj). 

The histograms in Figure [4] demonstrate that Main 
Pulses occasionally occur well above 5 GHz, even though 
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Fig. 4.— Phase-resolved logarithmic histograms of the occurrence of single pulses. The histograms have been aligned so that the most 
probable occurrence of the Main Pulse is at pulse phase zero (0°). The logarith m of the n umb er of occurrences is normalized to the 
maximum for plotting. Data sources: v < 5 GHz, recorded at the Arecibo Telescope |Cordes et al.| (|2004|); v > 5 GHz, this work, recorded 
at the GBT. 


they are too infrequent to be detected in the intensity- 
weighted mean profiles in Figure [Tj In particular, Main 
Pulses are evident at 9 and 14 GHz, and a few were de¬ 
tected at 28 GHz. Interpulses occur up to 43 GHz, the 
highest frequency we used. In Hankins & Eilek (20071 
are shown individual examples of Mam Pulses captured 
with our UHTRS at 9 GHz; in Paper 2 we show examples 
of single Main Pulses and Interpulses detected above 10 
GHz, including one Main Pulse recorded at 43 GHz. We 
note that both Main Pulses and Interpulses are so rare 
and weak at 43 GHz that the mean profile does not show 
either component; thus we did not include 43-GHz data 
in Figure |T| 

The histograms in Figure [4] show no significant detec¬ 
tions at the phases of any of the other pulse compo¬ 
nents. The non-zero histogram values outside the nomi¬ 
nal phases of the Main Pulse and Interpulses result from 
noise spikes or interference that exceeded the signal-to- 
noise ratio threshold. 


4.2. Single Pulse Occurrence Rates 


The pulse-count histograms in Figure [4] clearly show 
the striking difference between the count rates of the 
Main Pulse and those of the two Interpulses. At low 
frequencies the Main Pulse dominates, but above ~ 5 
GHz—the frequency range where the Interpulse changes 
nature—the Interpulse dominates the count rates, just 
as it does the intensity-weighted mean profile in Figure 

ffl 

This can also be shown in terms of count rates, ex¬ 
pressed as pulses/second. In Figure [5] we combine lower- 
frequency results from the literature with our current 
high-frequency data, to show the measured rates of the 
Main Pulse and two Interpulses, Amp and iVip. Al¬ 
though the measured count rates depend strongly on 
the observational systematics, such as telescope sensi¬ 
tivity and the Crab Nebula background level, the ratio 
Aip/Amp is much less sensitive to these effects. 

Figure [5] shows that bright Main Pulses are somewhat 
more frequent than bright Low-Frequency Interpulses be¬ 
low ~ 5 GHz, but the arrival rates differ only by a factor 
of a few. However, above ~ 5 GHz, the story changes. 
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Fig. 5. — Top: single pulse occurrence rates for Main Pulses, 
IVmp, and for Interpulses, A)p . Bottom: the rate ratio Nip /Nmp , 
comparing the occurrence rate of Interpulses to the occurrence rate 
of Main Pulses. Below 5 GHz Nip describes the occurrence rate 
of Low-Frequency Interpulses; above 5 GHz it describes the occur- 
rence rate of High-Frequen cy Interpulse s. Sources: 0.1 to 0.2 GHz , 
|Karuppusamy ct al. (201 21; 1.197 GH z.lPopov &; Stappersl (12007b : 
1.38 GHzl KaruppusamyTt al. | (|2010 I; 0.43, 1.18, 1.475 to 8.8 GHz, 
|Cordes et al. | ( |^!UU4[ ; 9, 14, 2U, 24, 28, 43 GHz, this work. 

Bright Main Pulses become increasingly rarer; we typi¬ 
cally caught only one every couple of minutes. However, 
High-Frequency Interpulses are much more frequent; we 
typically detected them 30-1000 times more often than 
Main Pulses in our continuously sampled observations 
between 9 and 43 GHz. This statistic, TVjp/IV mp, is some¬ 
what sensitive to the amplitude distr ibut ions of the sin¬ 
gle pulses (discussed below in Section 4.4 ); but when the 
rare Main Pulses above 10 GHz are detected, they are 
certainly no weaker than the High-Frequency Interpulses 
and sometimes much stronger. 

4.3. Durations of Single Pulses 

Pulse-count histograms (Figure El) and mean profiles 
(Figure]!]) do not tell the full story ofthe emission physics 
in the Crab pulsar, because individual pulses emitted at 
the phase of a given component are generally much nar¬ 
rower than the ensemble width of the component. Our 
UHTRS observations, obtained at time resolutions from 
100 /is down to 0.25 ns, allow us to resolve individual 
bright pulses, and thus characterize their duration, over 
a wide frequency range. We have recorded single pulses 
at the VLA at 0.33, 1.4, 4.8 and 8.4 GHz with 10- to 
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Fig. 6.— The averaged equivalent widths of the pulse intensi¬ 
ties at various frequencies are shown. The Main Pulse equivalent 
widths are fitted with a power-law, EWmp <x v~ 1-76±0.13^ between 
1.24 and 9 GHz. At each of the frequencies 14, 20, 28 and 43 GHz 
only one Main Pulse was captured by the UHTRS. For these widths 
no uncertainty bars are shown. At 24 GHz two Main Pulses were 
detected. The Low-Frequency Interpulse fit is, EWlfip °c v~ 1 ■ 38 , 
between 1.4 and 5 GHz. Because this fit has only two degrees of 
freedom, no error estimate is given. The High-Frequency Interpulse 
fit is EWhfip oc i/ _1 -18±0.06 between 6.5 and 24.25 GHz. No sin¬ 
gle High-Frequency Interpulses were detected by the UHTRS above 
28 GHz. For comparison with widths of the mean-profile compo¬ 
nents, note that 1 /rs ~ 0.1° of phase. 

intensity autocorrelation function (ACF) give similar re¬ 
sults, differing only by a small multiplicative factor. (The 
intensity ACF shows a “zero-lag spike” due to the noisy 
structure of the pulse and the contribution of the receiver 
and sky noise. Therefore, the ACF zero-lag amplitude is 
strongly dependent upon the signal-to-noise ratio (Rick- 
ett 19751. If used in the calculation of ACF equivalent 
width, tne zero-lag value would strongly bias the width 
estimate to narrower values. Therefore, for our analysis 
we truncated the ACF value at zero lag and replaced it 
by fitting a parabola to adjacent lags before computing 
the ACF equivalent width.) In Table El we present the 
number of pulses captured with the UHTRS at each fre¬ 
quency and the mean intensity equivalent width at each 
frequency. In Figure [6] we show the the equivalent width 
of the pulse intensity against frequency, averaged over 
all pulses recorded at each frequency. Below 2 GHz our 
mean equivalent widths for the Main Pulse are consis- 
tent with the range of si ngle pulse widths p resented in 
Crossley et al. (2010) and |Majid et al. (2011). 


ley et al. 2010) 

; at the Arecibo telescope with fret 

luency 

profile components (Table |3] and Figure |) 

bands 4-6, 6-8, 

and 8-10.5 GHz and 0.4-ns time resolu- 

result: individual bright pulses are rnucl 


tion (partial results reported in Hankins & Eilek 2007); 
and at the GBT with frequency bands 8-10, 12-15.6, 18- 
22.4, 22-26.5, 28-32 and 41-46 GHz, with 0.25-ns time 
resolution (as reported in this paper, and in more detail 
in Paper 2 to follow). 

Because the time signature of individual Main Pulses 
is quite diffe rent from that of indivi dual High-Frequency 
Interpulses (|Hankins & Eilek 2007), we need a measure 
of width that is not overly sensitive to the detailed pulse 
structure. We choose to calculate the equivalent width 
(EW) of the pulse intensity. After experimentation, we 
found that the EW of the intensity and the EW of the 


the duration of th eir associated mean-profile component. 
From Section 13.31 we know the width of the Gaussian fit 
to the mean profile of the Main Pulse, at a few GHz, 
is ~ 370 \js. For comparison, single bright Main Pulses 
last at most only ~ 20 /xs, and usually less, over the 
same frequency range. Similarly, the width of the High- 
Frequency Interpulse in the mean profile is ~ 730 /xs 
between 5 and 30 GHz, but the width of individual bright 
Interpulses is no more than ~ 2 p,s, and usually less, in 
this frequency range. This result makes it clear that 
the mean-profile components do not arise from spatially 
extended regions with uniform radio emissivity. Rather, 
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TABLE 4 

Single Pulse Equivalent Widths 


Frequency 

(GHz) 

Main Pulse 
Width (/is) 

Main 

Pulses 

Interpulse 
Width (/is) 

Inter¬ 

pulses 

0.333 

409. ± 51. 

47 

343. 

± 7. 

29 

1.241 

23. ± 

12. 

32 




1.385 

11. ± 

9. 

39 




1.435 

28. ± 

9. 

386 

11. 

± 9. 

35 

1.714 

10. ± 

7. 

76 




3.465 

1.8 ± 

1.4 

11 




4.150 

2.1 ± 

2.0 

27 




4.500 

1.5 ± 

1.1 

15 




5.013 

1.8 ± 

1.6 

79 

2.0 

± 0.3 

3 

5.500 

1.6 ± 

1.4 

45 




6.500 

1.05 ± 

0.9 

43 

1.7 

± 0.5 

33 

7.300 

0.74 ± 

0.5 

21 

1.4 

± 0.5 

40 

9.250 

0.70 ± 

0.5 

88 

1.3 

± 0.4 

190 

14.000 

0.32 


1 

0.64 

± 0.2 

54 

20.000 

0.53 


1 

0.41 ± 0.1 

130 

24.250 

0.65 


2 

0.38 ±0.1 

43 

28.000 

0.82 


1 

0.30 ± 0.07 

16 

43.250 

0.77 


1 





the mean profile components are an envelope inside of 
which short-lived bursts occur. 

Figure [6] also shows that the bursts which constitute 
single pulses become shorter at high frequencies. This 
is true for Main Pulses, Low-Frequency Interpulses and 
High-Frequency Interpulses. To quantify this behavior, 
we carried out least-squares fits of the form EW(^) oc v a . 
For Main Pulses between 1.24 and 9 GHz, we find a = 
— 1.76 ± 0.13. This fit is app roximately consistent with 
the value of a ~ —2 given by Crossley et al. (2010) over 
a small frequency range (see their Figure 6). For Low- 
Frequency Interpulses in the range 1.4 to 4.8 GHz we 
find a = —1.38. For High-Frequency Interpulses in the 
range 6.5 to 24.25 GHz we find a = —1.18 ± 0.06. This 
fit is approximately consistent with the widths reported 
for sin gle High-Frequency Intcrpulses by |Hankins fc Eilek| 
(2007) over a much narrower frequency range. 

We also computed the intensity equivalent widths of 
Main Pulses we captured at 0.33 GHz, finding the aver¬ 
age width to be 400 ± 50 ps. At this frequency the pulse 
shape is dominat ed b y interst ellar and Nebular scattering 
broadening (e.g., Crossley et al. 2010). We can, however, 
estimate the intrinsic single-pulse width at this frequency 
by extrapolating our power-law fit (upper panel of Figure 
[6| down to 0.33 GHz, which gives ~ 200 ps. Scaling the 
typical sc atteri n g broadening time of tq(1 GHz) ~ 2 ps 
at 1 GHz (Crossley et al.|2010 Figure 6) by v~ A down to 
0.33 GHz we obtain td(0.33~GHz) ~ 170 ps. If the 200- 
ps intrinsic pulse is convolved with a 170-/rs scattering 
function then the expected observed pulse width should 
be on the order of \/200 2 + 170 2 ~ 260 ps, which is not 
inconsistent with our width measurement. 


4.4. Fluence Distributions 

We used the continuously sampled GUPPI data to de¬ 
termine the fluence distribution of High-Frequency In¬ 
terpulses at each frequency we observe between 9 and 
43 GHz, averaging over 2-second time bins. (There were 
not enough Main Pulses in our data to provide reliable 
statistics). As is conventional, we fitted our fluence dis¬ 
tributions with power laws. In Figure [TJ also in Table 
[5j we show the distribution of power-law indices of our 
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Fig. 7. — The power-law indices of single-pulse fluence distribu¬ 
tions. The data for is > 9 GHz are for High-Frequency Interpulses 
measured in this work; other points have been extracted from the 
literature. Below 5 GHz the power-law index describes the flu¬ 
ence distribution of Low-Frequency Interpulses; above 5 GHz it 
describes the fluence distribution of High-Frequency Interpulses. 
Sources for the data points are given in Table [5] Uncertainty error 
bars are shown when available. 


TABLE 5 

Power-Law Indices of Fluence Distributions 


Frequency 

(GHz) 

Main Pulse 
Index 

Interpulse 

Index 

Source 

0.117 

-2.12 ±0.09 

-2.70 ±0.29 

(a) 

0.139 

-1.65 ±0.14 

-1.20 ±0.13 

(a) 

0.142 

-1.51 ± 0.05 

-1.60 ±0.09 

(a) 

0.148 

-1.78 ±0.10 

-1.60 ±0.12 

(a) 

0.157 

-1.83 ±0.10 

-0.70 ±0.05 

(a) 

0.163 

-2.38 ±0.15 


(a) 

0.174 

-2.39 ±0.12 

-2.05 ± 0.16 

(a) 

0.146 

-2.50 


(b) 

0.330 

-2.50 ±0.2 

-2.40 ±0.2 

(b) 

0.330 

-2.95 ±0.09 

-3.1 ±0.2 

(c) 

0.430 

-2.3 


(d) 

0.812 

-3.46 ±0.04 


(e) 

1.200 

-2.1 ±0.3 

-2.4 ±0.4 

(c) 

1.200 

-3.1 ±0.2 

-2.81 ± 0.03 

(c) 

1.373 

-2.79 ±0.01 

-3.06 ±0.06 

(f) 

1.435 

-3.06 ±0.00 

-1.70 

(g) 

8.800 


-2.9 

(d) 

9.375 


-4.00 

(h) 

14.375 


-2.65 

(It) 

20.375 


-3.64 

(h) 

24.375 


-4.65 

(h) 

24.874 


-4.07 

(h) 

28.375 


-5.95 

(h) 

28.375 


-4.38 

(h) 

28.375 


-3.77 

(h) 


N ote. — Sources: ( a) Karuppusamy e t al.l J20121); 
(b) l A rgyle fe Gower (1972); (c) Mickaliger et al. 
i 2012): (d) Uordes et al. (2U04); (e) Lundgrcn et al. 
i 1995): (f) Karuppusamy et al. (2012); (g) Moltett & 
ifaBsI pSS) ; (h) this work. - f '- 


data from 9 to 30 GHz along with power-law indices from 
lower frequency observations from the literature. 

There are, of course, many effects that can bias the in¬ 
dices. In addition to different systematic effects which 
can change with different observations, and the diffi¬ 
culty of establishing true power-law behavior when only 
a small range of fluence is sampled, the fluence distri- 
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Fig. 8.— Occurrence statistics of the Main Pulses and High-Frequency Interpulses for 9 < v < 28 GHz. The left-hand plots show the 
2-second average fluences as a function of time. The central panels show the correlation between fluence and pulse occurrence rates with 
the correlation coefficient. The right-hand panels show histograms of the fraction of the 2-second intervals during which a particular pulse 
rate occurred. The rate histogram ordinate exceeds unity for the intervals in which no pulses were detected. 


bution can vary with frequency and between the Mai n 
Pulse and the two Interpulses. For example, |Popov fc 


Stappers (2007), working at 1.197 GHz, found that when 


the giant pulses are sorted according to their width, the 
Main Pulse distributions are best fit by two power- law 
indices, and the Interpulse by one. Also Karuppusamy 


et al. (20101, working at 1.38 GHz, found the amplitude 
distributions are power-law only at the higher intensity 
tail of the distributions. Therefore, we do not interpret 
these power-law indices as proof of true power-law be¬ 
havior of the fluence distribution, but rather as a general 
indication of the relative abundance of stronger or weaker 
pulses in that distribution. 

Nonetheless, despite these uncertainties, Figure [7] 
shows that the Main Pulse and b oth Interpulses share 

([20()4|) 
ativeiy 


a m odest trend. As also n oted by Cordes et al. 
and Crossley et al. (20101, there seem to be re. 
fewer “super-giant" pulses at higher frequencies in each 
of these components. 


4.5. Bursts of Bright Pulses 

We also used the continuously sampled GUPPI data 
to study the arrival statistics of bright single pulses, 
at each observing frequency between 9 and 28 GHz. 
The different frequencies were not simultaneous, because 


they were observed on different days, but each Main 
Pulse/Interpulse time series at a given frequency was ex¬ 
tracted from data taken on one observing run. In Fig¬ 
ure [8] we present three metrics to characterize the arrival 
statistics on timescales from 2 to 10,000 seconds. 

4.5.1. The data: Fluctuating Arrival Rates and Fluences 

In the left panel of Figure [8] we show the time series of 
measured fluences at the phases of the Main Pulse and 
High-Frequency Interpulse, averaged into 2-second bins. 
This figure shows that single pulses at these frequencies 
do not arrive steadily, but rather occur in bursts last¬ 
ing from a few seconds to several minutes. The sepa¬ 
ration between bursts is highly variable; sometimes the 
bursts are separated by only ~ 100 seconds, but often 
it is much longer. There have been observing sessions 
where we have seen no single pulses for several hours. 
This behavior differs from previous observations at lower 
frequencies, where pulses are virtually always detectable 
in any one-second interval (for instance, see Figure [5]). 

The middle panel of Figure [8] presents the correlation 
between the pulse occurrence rates and the pulse fluence. 
A clear correlation is seen, demonstrating that brighter 
bursts tend to bunch together. The right panel of Figure 
[8] presents histograms of the pulse arrival rates, charac¬ 
terized as the fraction of the 2-second intervals in which 
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a particular pulse rate was seen. These histograms show 
that most of the time we recorded none or a only very 
few pulses per second, but that occasional multi-pulse 
bursts were possible, especially for the High-Frequency 
Interpulse which dominates at these high frequencies. 

4.5.2. The Fluctuations: Intrinsic or Interstellar? 

Rapid fluctuations such as seen in the left panel of Fig- 
ure[8]are often ascribed to diffractive interstellar scintilla- 
tion (PISS) in th e strong scattering regime (e.g., Cordes 
& Rickett 19981. Because we have the intensity time 


series over a wide range of wavelengths, we can deter¬ 
mine whether DISS is the cause of the burstiness we see 
between 9 and 28 GHz in the Crab pulsar. 

In Figure[9]we show the ACFs of each time series shown 
in the left panel of Figure [8] The DISS timescale, tdiss, 
is ty pically defin ed as the half-power point of this ACF 
(e.g., Cor dgs et al. ;;2004). Direct inspection of the ACFs 
in Figure [SJshows that tdiss fluctuates around ~ 50 to 
100 s over this range, but does not show any system¬ 
atic increase with observing frequency v. This behavior 
is inconsistent with standard DISS theory (e.g., Rickett 
1990), which predicts that tdiss M °c v x , where x = 1.0 
for Gaussian turbulence. Thus, tdiss should increase by 
a factor of ~ 3 over our 9 to 28 GHz observed frequency 
range — which is not the case. This suggests that the 
burstiness we observe is intrinsic to the pulsar. 

We can also ma ke some simple scaling e stimates, fol¬ 
lowing methods in Cordes & Rickett (1998) also B. Rick¬ 
ett, private communication (2014). The exponential de¬ 
cay constant at 330 MHz, td ~ 200 qs , is th oug ht to 
be due to interstellar scintillation (e.g., |Crossley et ah 
2010). If this is the case, we can estimate the interstellar 
scintillation bandwidth at 330 MHz as Ai-'iss (330 MHz) = 
1/27 ttd ~ 800 Hz. This quantity is predicted to increase 
with frequency, as Suiss yS , where s = 2(x + 2 ), thus 
s = 4.0 for Gaussian turbulence. For instance, we ex¬ 
pect 5v diss ~ 0.34 MHz at 1.5 GHz, 670 MHz at 10 
GHz, and 55 GHz at an observing frequency of 30 GHz. 
At observing frequencies where our receiver bandwidth, 
BW, is comparable to or less than the scintillation band¬ 
width (e. g., fo'Diss > 0.2 (BW), Cordes & Rickett 1998), 
burstiness such as we observe could be due to interstellar 
scintillation. For our 800-MHz GUPPI bandwidth, fre¬ 
quencies >,7 GHz could be affected by scintillation. How¬ 
ever, in addition to the fact that our observed tdiss does 
not show the expected linear increase with frequency, our 
scaling estimates also suggest the observed burstiness is 
intrinsic, as follows. 

The correlation bandwidth allows us to predict the 
fluctuation timescale expected from DISS, as tdiss iy) — 
(rp /v)(Sviss/v ) 1 / 2 , where tf = {cz/Zitv) 1 / 2 is the Fres¬ 
nel scale, z is the distance to the scattering screen, and 
v is the speed of the pul sar relative to that screen (e.g., 
Cordes & Rickett||1998|). We estimate v ~ 120 km/s, 
from an HST meas urement of the pulsar proper motion 
(|Kaplan et al.|2007), predicting tdiss ~ 27 s at 330 MHz, 
and ~ 80 s at 1 GHz. For a sanity check, we can compare 
our ex pected fl uctuation timescale at 1.5 GHz to the re¬ 
sults of Cordes et al. (2004]). These authors formed a time 
series ol giant pulses detected from the Crab pulsar, and 
measured the 1/e width of its ACF as ~ 25 s. Identify¬ 
ing this timescale as tdiss (which they argued should be 
detectable at 1.5 GHz with their observing bandwidth), 


they noted their result is a factor ~ 3 shorter than ex¬ 
pected from interstellar scintillation. They suggested the 
discrepancy is due to unusually high speeds of filaments 
within the Crab nebula; but intrinsic variability from the 
pulsar on shorter timescales than tdiss, may be another 
possibility. 

With these scaling arguments, we predict tdiss ~ 800 
s at 10 GHz, and ~ 2500 s at 30 GHz. Because these 
timescales are significantly longer than the burstiness we 
observed, we believe the burstiness we observe is not due 
to interstellar scintillation but rather is intrinsic to the 
radio emission mechanism in the Crab pulsar. 
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Fig. 9.— The normalized autocorrelation of the occurrence rates 
is shown. The rate values are the number of pulses per two-second 
interval. Where the ACF crosses the dashed line at 0.5 gives a 
characteristic time scale of the burstiness. 


5. SUMMARY AND DISCUSSION 

In this paper we have presented new, high-time- 
resolution observations of the Crab pulsar, at frequencies 
between 9 and 43 GHz. We combined our new work with 
previous data taken by our group and from the literature 
to characterize the mean radio profile and single-pulse 
statistics of radio emission from the pulsar at these high 
frequencies. 

5.1. The Split Personality of the Crab Pulsar 

Our main result is the dramatic change of the radio 
emission around frequencies of a few GHz, as follows. 

5.1.1. Low-frequency Radio Profile 

Below ~ 5 GHz, the Crab’s mean profile is dominated 
by two bright components, the Main Pulse and the (Low- 
Frequency) Interpulse. (Two weaker radio components, 
the Precursor and the Low-Frequency Component can 
also be detected below ~ 1 GHz, close to and leading 
the Main Pulse in phase.) The Main Pulse and Low- 
Frequency Interpulse are approximately phase-coincident 
with the two strong peaks in the pulsar’s high-energy pro¬ 
files. This near-coincidence of emission phases suggests 
the radio and high-energy emission occur in the same 
regions of the star’s magnetosphere, but that the radio 
emission sites are more spatially localized (or more tighly 
beamed) than the high-energy emission regions. 
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5.1.2. High-frequency Radio Profile 


The mean profile is dramatically different above ~ 5 
GHz; one would hardly believe one is looking at the same 
star. The Main Pulse disappears almost entirely above 
~ 5 GHz. The Interpulse — which is now the High- 
Frequency Interpulse — dominates the profile. It appears 
at a slightly earlier rotation phase than its low-frequency 
counterpart, but it still sits within the broad profile of 
the y-ray Interpulse. The two High-Frequency Compo¬ 
nents are also strong in the radio mean profile from ~ 4 
to 28 GHz. There is no sign of any component in op¬ 
tical, X-ray or 7 -ray mean profiles at the phases of the 
two High-Frequency Components, except for a possible 
(2.3<r) detection of the second High-Fre quency Compo - 
nent in the ,> 10 GeV mean profile of Abdo et al. (2010). 


5.1.3. Constraints on Future Models 

We do not, unfortunately, have a ready explanation for 
the geometrical or physical origin of the multiple radio 
components we see in the Crab pulsar. We therefore 
close by summarizing the key points about the different 
mean-profile components in the Crab pulsar, which we 
hope future models will address. 


5.1.4. How Sudden is the Change? 

The transition between low and high radio frequencies 
is sudden, but not “binary”. Although the Main Pulse 
disappears from the mean profile above 5 GHz, single 
bright Main Pulses can still be detected — and recorded 
at sub-ns time resolution up to 43 GHz. The change 
in the Interpulse may be sharper. We have not captured 
any single Low-Frequency Interpulses above 5 GHz, nor 
any High-Frequency Interpulses below 4 GHz. We have, 
however, seen single examples of both Interpulses — sep¬ 
arated by a few degrees of phase — in UHTRS dat a taken 
around 4 GHz. As we suggested in Section |3.2| the in¬ 
termediate phase of the mean-profile Interpulse at 4.15 
GHz may reflect a mixture of both types of Interpulses 
at that frequency. 


Frequency Interpulse by 7.3° ~ 640 /xs, and so it leads 
the maximum of the high-energy Interpulse by ~ 310 /its. 

5.1.6. Unsteady Radio Emission 

Radio emission from the Main Pulse and the High- 
Frequency Interpulse is quite unsteady. Our continuously 
sampled data show characteristic burst times ~ 30 to 50 
s at all frequencies we recorded between 9 and 28 GHz. 
We also found that bright bursts can bunch together; 
and, conversely, the star can go “quiet” (no strong radio 
bursts) for at least several hours. While we cannot def¬ 
initely prove that the burstiness is not interstellar, the 
frequency independence of the burst times suggests we 
are seeing fluctuations intrinsic to the emission process. 

5.1.7. Mean-profile Components Contain Short-lived Bursts 

Components of the mean radio profile are built up from 
the sum of many short-lived radio pulses. Our UHTRS 
data show that bright single pulses at the phases of the 
Main Pulse and both Interpulses typically last no more 
than a few microseconds. This is much shorter than the 
duration of the mean-profile components, which last sev¬ 
eral hundreds of microseconds. In addition, we show in 
Paper 2 that more than one bright single Interpulse can 
occur during one stellar rotation, within the phase win¬ 
dow of the mean-profile component. Thus, a mean-profile 
component does not represent steady radio emission dis¬ 
tributed over an extended region in the magnetosphere; 
rather, it is the envelope of where individual pulses can 
occur. 


5.1.8. Single Pulses at High Time Resolution 

Much more can be learned by studying single pulses 
at sub-ns time resolution. In Hankins & Eilek| (2007) 
we showed that the High-Frequency Interpulse has dif- 
ferent spectral, temporal and dispersion characteristics 
from the Main Pulse. This suggests that different phys¬ 
ical conditions exist in the emission regions, and/or the 
propagation paths, for each component. In Paper 2 we 
extend our single pulse analysis to higher radio frequen¬ 
cies and additional components — and demonstrate that 
the radio emission physics from the Crab pulsar is even 
more complex than we have seen up to now. 


5.1.5. Relation of Radio and High-energy Components 


The Main Pulse and both Interpulses occur at approx¬ 
imately the same phases as their counterparts are seen at 
optical, X-ray and 7 -rays, which suggests they come from 
similar regions in the star’s magnetosphere. There are, 
however, important differences. Both the Main Pulse 
and the Low-Frequency Interpulse are significantly nar¬ 
rower in phase than their high-energy components, and 
each of these radio components lags the corre spond ing 
high-energy maximum by ~ 20 0 - 300 /is (e.g., Abdo et 


al. 


2010 Zampiere et al. 2014| and references therein). 


Different timing results from different authors and tele¬ 


scopes display some scatter, but all agree that the radio 
pulses lag the high-energy pulses by a fraction of a ms. 
The High-Frequency Interpulse also sits within the broad 
peak of the high-energy Interpulse, and is also narrower 
than its high-energy counterpart, but it leads the Low- 
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